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A  GLOBAL  EXISTENCE  AND  UNIQUENESS  THEOREM  FOR  A  RICCATI  EQUATION* 

J.  COLBY  KEGLEY 


Abstract.  It  is  proved  that  a  Riccati  differential  equation  of  a  particular 
form  has  a  unique  solution  satisfying  the  conditions  that  it  is  to  exist  for 
large  values  of  the  independent  variable  _t_  and  to  have  its  graph  stay  above  a 
certain  line  for  large  t  .  It  is  then  proved  that  the  solution  exists  for  all  t  . 

■X/  Tf 

Two  forms  of  the  solution  are  developed  in  terms  of  the  confluent  hypergeometric 
functions.  An  application  of  these  results  is  made  to  an  asymptotic  stochastic 
analysis  of  a  noisy  duel  problem. 
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1*  Introduction.  This  paper  investigates  a  particular  form  (1)  in  section 
2  of  a  Riccati  equation  that  is  quadratic  in  the  independent  variable  t.  The 
approach  to  the  problem  of  the  existence  of  a  global  solution  is  not  from  the 
usual  initial-value  standpoint,  but  is  based  on  a  desired  feature  of  the  solution 
for  large  t  which  is  given  by  properties  (i)  and  (ii)  of  the  Theorem  in  section 
2.  The  form  (1)  of  the  equation  makes  it  easy  to  draw  rough  sketches  of  how  the 
solutions  behave  depending  on  where  the  initial  point  is  selected.  In  particular, 
it  becomes  plausible  that  there  exists  a  solution  defined  for  all  _t_  that 
satisfies  properties  (i)  and  (ii),  but  it  is  by  no  means  clear  that  there  is  only 
one  such  solution.  That  this  is  the  case  indicates  that  this  distinguished 
solution  is  extremely  unstable.  Indeed,  one  of  the  implications  of  Lemma  5  in 
section  3  is  that  every  other  solution  diverges  from  the  distinguished  solution 
as  t  -►  00 . 

Our  investigation  is  motivated  by  the  approach  used  in  [3]  and  [6]  to 
analyze  the  equal-accuracy  noisy  duel  problem  for  two  players  having  finite 
unequal  units  of  ammunition.  This  approach  leads  to  asymptotic  distributions 
of  normalized  times  of  first  fire  for  the  two  players.  The  hazard  rates  for 
these  distributions  are  expressed  in  terms  of  a  solution  to  a  Riccati  equation 
of  the  form  (1),  and  the  distributions  themselves  are  expressed  in  terms  of  a 
solution  to  a  related  Hermite  equation. 

A  brief  outline  of  these  connections  is  given  in  section  4.  The  reader 
may  find  it  helpful  to  read  that  section  in  conjunction  with  the  statement  of 
the  Theorem  to  understand  the  reason  for  deriving  the  various  properties  of  the 
distinguished  solution. 


2.  Statement  of  the  Theorem.  The  principal  conclusions  we  desire  can  be 


stated  as  follows. 
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THEOREM.  Suppose  a  is  a  positive  number  and  <t>2  are  linear  functions 

4>^(0  =  where  <  6^  and  g^  >  0.  Then  there  is  exactly  one  solution 

of  the  Riccati  equation 


(!)  v'(t)  =  a[v(t)  -  2(J,1(t)][v(t)  -  24>2(t) ] 

t_hat  has  the  following  two  properties:  There  is  a  number  t q 
(i)  The  domain  of  v  includes  the  interval  [tg,°°); 
(ii)  v(t;  201(t)  >  0  for  t  _>  tQ. 

Moreover,  thi ■  •  I ution  has  the  additional  properties : 

(iii)  The  domain  of  v  _is  (-00,00). 

(iv)  v ( t )  -  2<J>1(t)  >  0  tor  aU  I* 

(v)  /°°[v(t)  -  2<J>1(t)]dt  «  /°°[v(t)  -  2<(»2(t)]dt  = 

i)  v(t)  -  2<t> ^ (t )  -*  0  as  t  -*■  <*>. 

'  ii)  v'(t)  -*■  2 0^  as  t  -*■ 

If  ri  addition,  B2  <  0,  then  the  following  hold: 
ii)  1^  tp  ^s  a  number,  then  the  conditions 

v (t )  -  2<t2 ( t)  >  0  for  t  >  tQ 
and  v(tQ)  -  2^2(t0)  -  0 
hold  exactly  when 

*^1  ( Cq)  “  ^2^0^  “  Z0l/^l“^2 ^2a’ 


such  that 
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where  z^  is^  the  real  zero  of  the  Weber  parabolic  cylinder  function 
Dp+1  with  P  “  B1/(B2-Bi). 

(ix)  With  t^  as  in  (viii)  arid,  for  i  -  1,2,  we  define  fj(t)  -  x1(t)$1(t), 
where  x^t)  -  n[v(t)  -  2^(0]  and  *t(t)  -  exp(-J^  x^Tjdi) ,  have: 
(ixa)  ( t )  <  $2(t)  when  tQ  <  t  <  -tQ  -  2n  and 
♦  ](t)  >  $2^^  when  t  >  -tQ  -  2n. 
where  t  ■  -n  the  solution  x:(t)  -  x2(t). 
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(ixb)  is  decreasing  and  positive  on  (-00  ,  00)  . 

(ixc)  is  positive  on  (t^,00)  and  has  a  maximum  value  that  occurs 

-a— number  >  -  n  >  tg. 

(ixd)  f  t f ~ (t)d t  <  00  ■ — - - 

(ixe)  C  tf ,  (t)dt  =  °° 
c0 

n  the  process  of  proving  this  theorem,  two  forms  of  the  solutions  are 
dev  loped . 


(A) 

where 

and 

with: 


V(t)  =  zm*0  ■ 


s  =  6(t+r|),  6  =  /otCB^-f^)  >  r)  =  (Yj~Y2) 

^(s)  =  t,  yQ(s)  +  y^s) , 

C  =  -2T (a  +  1/2)  /r (a)  ,  a  = 


Yq(s)  =  s  +  l/2>  3/2;  s  ), 

yj(s)  =  1/2;  s2) , 


and  denotes  the  confluent  hypergeometric  function 

p  /  h  ^  -  T(b)  r  r(a-Ki)  ^ 
1F1 (a ’  b’  z)  "  F(a)  n^0  T(b+n)  n!  * 


(B) 


where 


v(t)  -  2«2(t)  . 

z  =  6/2(t+ri),  p  *  Bj/(B2*'$p, 


and  Dp  is  the  Weber  parabolic  cylinder  function  (cf.  [7l)» 

D  (z)  ■  r(l/2)2P^2exp(-z2/4)R  (z)  , 

P  P 

wltb  Rp(z)  ,F,C-p/2,  1/2;  z2/2)  -  f^2)  ih <> «  ‘  •  3/2> 


■gszamam 
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3.  Proof  of  the  Theorem.  The  demonstration  of  the  conclusions  is  broken  down 
into  several  stages. 

LEMMA  1.  A  function  v  is  a  solution  of  equation  (1)  on  an  interval  I 
exactly  when 

_  1  x* 

(2)  "  1~ 

where  x(t)  ^  0  for  t  in  I  and  is  a  solution  of 

(3)  x"(t)  -  q (t)x(t )  =  0 

2  2 

with  q  =  a  + 

Proof.  We  re-write  equation  (1)  in  the  form 

v'  +  2Av  +  Bv2  -  C  =  0 

where  A  =  a(<fr  ,  B  =  -a,  and  C  » 

We  then  apply  the  result  in  Reid  [5]  that  v  is  a  solution  of  (1)  on  an  interval 
I  if,  and  only  if,  v  *  u/x,  where  x(t)  £  0  on  1^  and  the  pair  (x,u)  is 
a  solution  on  I  of  the  linear  system 

x*  *  Ax  +  Bu 

(4) 

u ’  =  Cx  -  Au. 

But  this  system  is  equivalent  to  equation  (3),  as  can  be  seen  through  the  connec¬ 
tion  u  *  (xf-Ax)/B.  Calculating  v  *  u/x  then  gives  the  form  (2). 

In  order  to  transform  equation  (3)  into  more  comprehensible  forms,  first 
we  make  a  change  of  independent  variable. 

LEMMA  2.  The  general  solution  of  equation  (1)  i£ 

»<t)  -  *l(t)  ♦  ♦jCO  -  i  . 

where  s  -  6(t+n),  <5  “  »  H  ■  •> 


and  w  is  a  non-vanishing  solution  of  the  Weber  equation 

(5)  w"(s)  +  (g-s2)w(s)  =  0 
with  €  =  (32+01)/(02-61)  • 

While  the  form  (5)  is  simpler  than  the  form  (3),  it  is  not  easy  to  see  when 

its  solutions  are  non-vanishing;  the  equation  has  arT os'c fix ax\~ri-y- y  a _1  ab out 

s  =  0  if  6  >  0.  However,  we  can  make  the  change  of  dependent  variable 
2 

y(s)  =  exp(s  /2)w(s),  which  transforms  (5)  into  a  Hermite  equation  and  clearly 
preserves  the  non-vanishing  of  solutions.  In  fact,  going  through  the  calculations 
gives  the  following  result. 

LEMMA  3.  The  general  solution  of  equat ion  (1)  rs 

(6)  v(t>-  2^,(0 

where  y  is  a  non-vanishing  solution  of  the  Hermite  equation 

(7)  yM(s)  -  2sy’(s)  -  4ay(s)  =  0 

with  a  =  31/2(B1-62) * 

Before  continuing,  we  point  out  that  the  procedure  of  transforming  an  equa¬ 
tion  of  the  form  (3),  where  q  is  quadratic,  first  into  the  form  (5)  and  then 
into  the  form  (7)  is  well-known.  It  is  used,  for  example,  in  solving  che  time- 
independent  Schrodinger  equation  for  a  harmonic  oscillator. 

Now,  the  general  solution  of  (7)  can  be  expressed  in  terms  of  the  confluent 
hypergeometric  functions.  In  fact,  we  have  the  following  result,  which  may  be 
verified  by  direct  calculation  or  by  referring  to  Slater  [7]. 

LEMMA  4.  Let  y^  and  y^  denote  the  solutions  of  equation  (7)  that  satisfy 
the  initial  conditions 
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y0(0)  -  o,  Yq(0)  -  1 

yx  (0)  =  1,  y|(0)  -  o. 

Then  the  functions  yQ  and  y1  are  given  by 

(8)  yQ(s)  =  s  lFl^a  +  ^2»  s2^’ 

_ 

( 9 )  "  y j-C&I.  =^F^(a, ^  1  / 2 ;  s2  ) . 

We  now  focus  our  attention  on  property  (ii)  in  the  statement  of  the  Theorem. 

The  next  result  shows  that,  up  to  a  multiplicative  constant,  there  is  only  one 
solution  of  (7)  which,  when  substituted  in  (6),  can  possibly  work. 

LEMMA  5.  Every  non-trivial  solution  y  -  cqYq  +  c^y^  2l  e9uat i°n  (?)  ^as  t he 
property  that  y'(s)/y(s)  «>  as  s  -*  qo  unless  the  constants  and  c^  satisfy 

the  relation 

cor(a)  +  2ClT(a  +  1/2)  =  0. 

Proof.  We  apply  two  results  about  confluent  hypergeometric  functions  given  in 
Slater  [7].  We  have  the  derivative  relation 
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y(s) 


r(l/2)exp(s2)s2a  1 
2T(a) T(a+l/2) 


{cor(a)[l40(s 


)] 


+  2Cir(a+l/2)[l+0(s  2)]|  . 


Therefore,  as  s  +  oo  we  have  y!(s)/2sy(s)  +  1,  which  implies  yf(s)/y(s)  ->  oo, 
unless  the  constants  and  c^  make  this  form  indeterminate.  But  this  is 

precisely  when  cQr(a)  +  2c^T(a+l/2)  =  0. 

Now  recall  that  the  variables  s  and  t  are  related  by  s  =  6(t4rj), 
where  6  >  0,  so  the  conditions  s  oo  and  t  -+  oo  are  equivalent.  Then  the 
form  (6)  for  v(t)  shows  that  a  solution  that  exists  for  large  t  will  have 
v(t)  -  2<J>2  ( t )  ->•  “°°  as  t  -*  oo  unless  c^  and  c^  satisfy  the  relation  stared 
in  Lemma  5.  Furthermore,  the  solution  of  (7)  enters  into  the  form  (6)  only 
through  the  ratio  yf/y,  so  one  of  the  constants  cQ  and  cj  may  be  chosen 
arbitrarily.  Since  we  want  to  avoid  solutions  y  that  vanish,  we  choose  c^  =3 
anc*  c0  =  -2F(a+L/2) /T(a) .  We  summarize  what  we  have  obtained  so  far  as  follows. 

LEMMA  6.  A  necessary  condition  for  a  solution  v  of  equation  (1)  to  ha'  ' 
properties  (i)  and  (ii)  of  the  Theorem  is  that 

(10)  v(t)  =  2<Mt)  -  -iirrT  > 

1  a  ip(s) 

where  4>(s)  =  C,yQ(s )  +yL(s)  with  q  =  -2r(a+l/2)/r(a)  and  y  ,y  are  defined 
by  formulas  (8)  and  (9) ,  respectively . 

Notice  that  (10)  is  the  form  (A)  of  the  solution  v_  that  is  given  in  the 
remarks  following  the  Theorem.  We  now  proceed  to  show  that  the  particular  solu¬ 
tion  ip  of  (7)  forces  the  corresponding  solution  v_  of  (1)  to  have  properties 
(i)  through  (ix)  of  the  Theorem. 

LEMMA  7.  The  solution  \p  =  QyQ  +  of  equation  (7)  satisfies  the  inequalities 
i|;(s)  >  0  and  (s)  <  0  for  all  s. 


i 


Proof.  The  major  theoretical  tool  we  need  to  prove  this  result  is  stated 

in  the  Appendix.  In  order  to  apply  that  theorem  to  our  problem,  we  put  equation 

2 

(7)  in  self-adjoint  form  by  multiplying  both  sides  by  exp(-s  ).  The  result  is 


the  equivalent  equation 


(ry * ) 1  -  py  =  0, 


where  r(s)  =  exp(-s  )  and  p(s)  =  4ar(s).  Since  r  and  p  are  continuous 
with  r(s)  >  0  and  p(s)  >  0  for  all  s_,  we  can  conclude  that  if  y  = 

is  the  solution  of  (7)  with  0  =  limg_^Q  -  yj(s)/yg(s),  then  y(s)  >  0  and 
y,(s)  <  0  for  all  js.  We  now  show  that  0  -  £. 

To  do  this,  we  proceed  as  in  the  proof*  of  Lemma  5.  Using  the  definition  of 
y^  and  y^  ,  and  the  asymptotic  expansion  of  the  confluent  hypergeometric  functions 
again,  we  obtain,  as  s  -*•  <»: 

yl(s)  -  ~Y(a)‘  eyp^s  U+0(s  )) 

and  yQ(s)  =  2F(a+f/2)exp^s2^s2a  lO+0(s  2))- 

This  makes  it  clear  that  0  =  -2r (a+1 /2) /V (a)  =  £  . 

By  referring  to  the  form  (10)  and  applying  the  result  of  Lemma  7,  we 
immediately  have; 

COROLLARY.  The  solution  v  of  equation  (1)  defined  by  (10)  sat isf ies 
properties  (iii)  and  (iv)  and,  a  fortiori,  satisfies  properties  (i)  and  (ii). 

In  order  to  tackle  properties  (v)  through  (viii),  we  develop  the  second  form 
(B)  of  the  solution  v. 


LEMMA  8.  The  solut ion  ip  *  of  equation  (7)  can  be  written  in  the 


,,  v  2t (a+i/2)  ,  2,,_  #  , 

*<■>  “  — r'(~l~2~)  exp(z  /4)Dp(z) 


9 


where  z  =  s/2  =  6/2(t+0)  and  is  the  Weber  parabolic  cylinder  function  with 

p  =  -2a  =  61/(B2“31). 

Proof.  The  result  follows  by  using  the  definition  of  and  simplifying 

the  right-hand  side  of  (11). 


LEMMA  9.  The  solut ion  v_  £f  equation  (1)  def ined  by  (10)  can  be  wri t  ten  in 
the  form 

d,+i <*> 

)  { 

P 

Proof.  To  obtain  this  form,  first  we  use  formula  (11)  for  ^  and  calculate 
Keeping  in  mind  that  z  =  s/2,  we  obtain 


(12) 


v(t>  -  *vo  +  ¥Hr^r- 


-  e  +  Jl  Dp(z)  . 

*(s)  dw 

P 


Then  we  use  the  Identity  (cf.  [A]) 


(13) 


D'(z)  =  (z/2)D  (z)  -  D  (z). 
p  P  p+1 


The  result  is  that 
(14) 


i'l?.)  =  2s  -  Ji  } 

<P(s)  D  (z) 


Substitution  of  this  expression  into  formula  (10)  and  use  of  the  relation 
s  =  6(t+n)  give  the  form  (12)  for  v,  which  is  the  form  (B)  that  was  claimed. 

Properties  (v) ,  (vi),  and  (vii)  can  now  be  attacked  by  using  the  following 
asymptotic  expansion  of  the  parabolic  cylinder  functions  (cf.  [4]).  As  z  <», 


(15)  Dp(z)  =  exP(“Z2M)  zP[l  -  +  0(z”^)]. 

2z 

LEMMA  10.  If  \p  =  ^y^  +  y^  ,  then  \p(s)  ■+  0*  as  s  ->  oo. 

Proof.  We  return  to  formula  (11)  for  keeping  in  mind  the  connection 
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z  *  s/T  and  the  fact  that  p  =  -2a  x  0.  Substitution  of  the  result  of  (15)  into 
(11)  yields,  as  s  -+  oo,  * 

(16)  4>(s)  =  2  7((t/2)2)  zP[1  ~  2S^r  + 

so  ijj(s)  ->  0  as  s  ->  oo.  Since  ijj(s)  >  0  for  all  s,  the  conclusion  follows. 


COROLLARY . 


The  solution  v  o£  equation  (1)  defined  (10)  satisfies  property  (v) , 


Proof.  Using  the  form  (10)  with  the  connection  s  =  6(t+n)  shows  that  if  we 
fix  some  t^  and  let  be  the  corresponding  value  of  s,  then 

(17)  f.  (v-2<(,  )  =  -  -J-  Is  (f/<J>)  -  “  ^*(s)/«s  )]. 

co  a  0  a 

sOO 

But  s  ->  oo  as  t  oo,  so  j  ( v—  2 tf> ^ )  =  00  follows  immediately  from  Lemma  10. 

to 

For  the  second  integral,  we  again  fix  t^  and  then  apply  property  (iv)  to  an 
interval  [tQ,t].  The  result  is 

fl  (v-2 <f>  )  >  ft  (2<J>,  ~2<p„) 

0  0 


=  ( ~e2  >  ( *■  2“to)  +  Oi“Y2)(t_to^ 


which  *>oo  as  t  since 


LEMMA  11.  If  ip  =  £y^  +  ,  then  t p'  (s)/ty(s)  0  as  s  ->  <». 

Proof.  We  begin  by  using  formula  (14),  recalling  again  that  z  =  s/2.  The 


result  is 


Jkl L*1  V2f*  -  -Vl(2)  i 

ii(s)  D  (z)  *' 

But  z  -*■  00  as  s  •+  ®>,  and  if  we  use  just  the  result 


Dp(z)  =  exp(-z2/4)  zP(l  +  0(z  2)) 


from  formula  (15),  we  obtain 


property  (vii). 
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Proof.  If  we  look  at  the  form  (12)  of  the  solution,  we  obtain 


*•<.)- *}<o  ♦*£(;*!)'<.>  £ 


Wfw- 


-  2B2  +  2  a- 


since  z  =  6/2(t+n).  But  z  00  as  t  *>  so  Lemma  12  implies  that 

(S2 

vf  (t)  2B2  +  2  —  as  t  -► 

Using  the  definition  6  =  then  gives  the  result. 

Next,  vje  use  the  following  result  about  the  parabolic  cylinder  functions. 
LEMMA  13.  If_  B2  <  0  <  3^  and  p  =  : 

(i)  D^(z)  >  0  for  all  z ; 

(ii)  D  ^  has  exactly  one  real  zero  zQ,  and  D  +1(z)  >  0  exactly  when 


z  >  zQ. 


Proof.  The  hypotheses  imply  that  0  <  p+1  <  1.  Hence,  the  result  follows 
immediately  (cf.  [1]). 

COROLLARY.  The  solution  v  of  equation  (1)  defined  by  (10)  satisfies 
property  (viii). 

Proof.  If  we  apply  Lemma  13  to  the  form  (12)  of  the  solution,  we  see  that 


v ( t )  -  2<J>2 (t)  >  0  for  t  >  tQ 
and  v(tQ)  -  2<|>2(t0)  =  0 

exactly  when  satisfies  zQ  =  6/2(t0+n)-  A  simple  calculation  using  the 

definitions  of  ^ ,  6  and  n  then  gives  the  result. 


J 
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LEMMA  14.  The  functions  def ined  in  (ix)  are  related  by 

(19)  ?2(t)  =  G (t)  ^(t), 

2  2 

where  G(t)  =  exp(s^  -  s  )  with  s  =  6(t+n)  and  =  SCt^-fn). 

2 

Proof.  Since  x^t)  =  cx[  v  (t)-2(}>^  (t )  ] ,  we  have  X£(t)  =  ^(t)  +  2($  (t+r0  > 
from  which  (19)  follows  easily. 

COROLLARY.  The  functions  satisfy  (ixa). 

Proof.  Since  t^  satisfies  (viii),  we  have 

<61_62) (t0+n)  =  W  '  MV  =  z0/(6rB2)/2a’ 

which  is  negative  because  the  zero  z^  of  with  0  <  p+1  <1  is 

negative  (cf.  [1]).  Hence,  tQ+r)  <  ®  since  6^  >  &2*  so  so  <  Therefore, 
G(t)  >  1  exactly  when  s^  <  s  <  -Sq,  i.e.,  when  t^  <  t  <  -t^  -  2n. 

LEMMA  15.  The  solution  \ jj  =  £y^  +  y^  £f  equation  (7)  has  l^"(s)  >  0 
for  all  s. 

Proof.  If  there  were  a  value  s^  at  which  \J/'(s^)  =  0,  then  by  taking 
the  derivative  with  respect  to  s  of  both  sides  of  (7)  with  y  replaced  by 
ip,  we  would  have  ^M,(s^)  *  (2+4a)^f  (s^) .  Since  just  the  hypotheses  <  B^ 
and  >  0  imply  a  >  0,  and  Lemma  7  implies  ^r(s^)  <  0,  it  follows  that 
(Sl)  <  0.  This  says  that  at  s^»  ip  changes  from  being  convex  to  concave, 
and  hence,  by  this  very  argument,  can  never  change  back  to  being  convex.  But 
that  contradicts  the  result  of  Lemma  7.  So,  ^"(s)  is  never  zero.  But,  by 
replacing  by  0  in  (7),  we  have  ^"(O)  =  4a  >  0.  Since  ip"  is  continuous, 
the  result  follows. 

COROLLARY.  The  function  f^  satisfies  (ixb). 

Proof.  Since  property  (iv)  implies  x^(t)  >  0  for  all  it  follows 

from  the  definition  of  f-  that  f^(t)  >  0  for  all  ^t.  To  show  that  f^  is 
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decreasing,  notice  that  form  (A)  of  the  solution  to  equation  (1) 
implies  that  x^t)  =  -Sip'  (s)  /ip(s)  =^riKs). 

Hence,  we  have 

(20)  ^(t)  =  Us)/UsQ). 

where,  as  in  Lemma  14,  s  =  6(t+n)  and  sQ  =  6(tQ+n).  But  then 

f^(t)  =  -5|(t)  =  -6ip'(  s)/(//(Sq),  so  f^(t)  =  -6^  ipn  (s)  /ip(s^)  <  0  by  Lemmas 

7  and  15. 

LEMMA  16.  The  function  satisfies  (ixc). 

Proof.  That  f2(t)  >  ^  f°r  t  >  £q  anc*  ^2^0^  =  ®  follows  from  the 

definition  of  f  2  and  (viii) .  Next,  we  show  that  f2(t)  ^  ®  as  t 

For,  the  definition  of  f^  and  relation  (19)  imply  f 2  *  -$2  =  “(G$^)f.  Using 

(20)  then  gives 

(2D  f2(t)  =  6G(t)*1(t)[2a-  ^(s)Ms)]. 

Since  t  >  t^  corresponds  to  s  >  s^,  (20)  and  Lemma  7  imply 

(22)  0  <  5  (t)  <  1  for  t  >  tQ. 

Also,  ^f(s)/ip(s)  ■>  0  as  t  -+■  00  by  Lemma  11.  Finally,  2sG(t)  0  as  t  « 

by  the  definition  of  Hence,  f9  is  a  continuous  function  with  f  2  ( t )  >  0 

on  (t0^  while  f2(tg)  =  ®  =  ^2^°°^  *  so  ^2  ^as  an  a^so^ute  maximum  in  (t^,00). 
To  facilitate  the  calculation  of  f^,  we  first  use  (14)  and  (19)  to  re-write 

(21)  as 

(23)  f2(t)  =  tS  $2(t)  Dp+i(z)/Dp(2). 

with  z  *  /26  (t+n)  as  before.  Taking  the  derivative  of  both  sides  of  (23) 
with  respect  to  ^  and  using  (18)  as  well  as  -  -f2>  we  find  that 

f£(t)»  2  62$2(t)[p+l  -  zDp+1(z)/Dp(z)]. 
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Since  t  >  corresponds  to  z  >  z^,  where  is  the  (negative)  zero  of 

Dp+^(z),  it  follows  from  Lemma  13  that  f^  (t)  >  0  for  z^  <  z  <_  0,  i.e., 

tg  <  t  _<  -n.  Hence,  the  maximum  of  occurs  at  some  >  -n« 

In  order  to  deal  with  (ixd)  and  (ixe),  we  use  f ^  and  integration 

by  parts  to  get 

(24)  JftTf,(T)dr  =  -  t$  (t)+  f*  $  (T)dx. 

0  0 

Then  (19)  and  (22)  easily  yield 

LEMMA  17.  The  function  f^  satisfies  (ixd). 

Finally,  (20)  and  the  asymptotic  expansion  expansion  (16)  for  ^  show 
that  $^(t)  behaves  like  tP  as  t  -*  where  p  -  6^ / ( &2_ )  *  Since 
^2  <  0  <  implies  that  -1  <  p  <  0,  it  follows  readily  from  (24)  that 

LEMMA  18.  The  function  f^  satisfies  (ixe). 

4.  An  application  to  a  noisy  duel  problem.  In  [3]  and  [6]  appears  a 
dynamic  programming  approach  to  the  vs.  equal-accuracy  noisy  duel 

problem,  where  the  positive  integers  m  <  n  represent  the  units  of  ammunition 
the  two  players  have.  The  approach  begins  by  allowing  either  of  the  players 
to  fire  a  unit  of  ammunition  only  at  times  corresponding  to  points  of  a  discrete 
grid  of  the  interval  [0,1],  which  is  interpreted  as  the  interval  of  probabilities 
of  either  player  destroying  the  other  if  a  unit  is  fired.  This  produces  a 
finite  sequence  of  simultaneous  games  whose  2*2  pay-off  matrices  are  deter¬ 
mined  by  proceeding  backwards  inductively  from  the  game  where  the  probability 
of  destruction  is  unity. 

Attention  is  focussed  on  an  interval  of  grid  points  at  which  the  players 
have  no  pure  strategy  and  which  surrounds  the  critical  probability  l/(m+n). 

It  is  found,  under  suitable  hypotheses  suggested  by  computer  implementation 
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of  the  above  approach,  that  the  value  of  the  game  in  this  interval  of  grid 
points  satisfies  a  difference  equation.  Dividing  both  sides  of  this  equation 
by  an  appropriate  normalization  factor  and  letting  the  mesh  of  the  grid  on 
[0,1]  approach  zero  leads  to  a  normalized  value  v  =  vm  n  the  ^>ame  that 

satisfies  a  Riccati  equation  of  the  form  (1)  on  the  interval  [-1,°°),  with 

2 

a  '  (m+n)  (nH-n-2) /(n-m)  , 

3.  =  m  c  2  / (m+n~ 1 )  ,  =  (-n/nOB.  , 


(m+n- 1 )  m,n 

~  2 (m+n)  c  . 

m  ,n-  i 


while  Yj  *  0  for  m  *  n  -  1, 


v  ,  (-1)  for  1  <  m  <  n- 1 , 
m ,  n- 1  — 


(m+n- 1 )  m ,n 


2  2 (m+n)  c 


m-1  ,n 


v  (-1)  for  1  <  m  <  n-1, 

ra-1 ,n  - 


while  ^2  *  0  for  m  =  1 , 


Here,  it  is  known  that  the  constants  c.  .  are  positive  for  i  <  j,  but 

cx  p  r  \o{ c  * 

analytic  expressions  for  these  constants  are  not  known^  However,  the  hypotheses 
ot  >  0  and  62  <  0  <  6^  are  evidently  satisfied.  Also,  it  is  established  in 
[3]  and  [6]  that  the  initial  condition  v(tg)  -  ^2^0^  =  ®  is  to  hold  when 
tg  *  -1.  But  this  does  not  seem  to  be  enough  information  to  attack  the 
existence  and  uniqueness  problem  for  (1)  on  [-I,00). 

Instead,  attention  is  turned  to  the  functions  defined  in  (ix)  with  tg  *  -1, 
which  corresponds  in  the  normalization  process  to  the  beginning  of  the  interval 
surrounding  probability  l/(m+n)  in  which  random  strategies  are  to  be  employed. 
The  functions  $^(t)  and  represent,  respectively,  the  probability 

that  the  weaker  player  and  the  stronger  player  has  a  normalized  time  of  first 
fire  occurring  at  or  after  The  functions  x^(t)  represent  the  corresponding 
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hazard  rates  for  the  cdf's  <K(t)  =  and  the  functions  f i ( t )  repre¬ 

sent  their  densities. 

One  of  the  facts  derived  in  [3]  and  [6]  is  that  the  weaker  player* s 
hazard  rate  x^(t)  is  to  be  positive  for  t  j>  -1.  Somewhat  surprisingly,  the 
assumption  that  the  solution  of  (1)  exists  for  large  and  that  x^(t)  be 

positive  for  large  produces  not  only  the  global  existence  and  uniqueness 

result  for  (1)  proved  herein,  but  also  some  properties  of  the  complements 
$^(t)  of  the  cdf’s  $^(t)  that  could  not  be  surmised  by  studying  the  computer 
runs  for  the  2  x  2  games,  namely: 

Property  (v)  implies  that 

J[_1xi(t)dt  =  °°  for  i  =  1,2, 

so  that  lim  $^(t)  =  exp  (~j_,  Xj  (t)dT  =  0, 

£-*O0 

which  implies  lim  $^(t)  =  1. 

t-+O0 

This  says  that  the  probability  is  unity  that  each  player  fires  at  some  time 
in  the  normalized  interval  -1  t  <  00  during  which  random  strategies  are 
employed , 

Property  (ixa)  states  that  the  probability  is  greater  not  only  for  the 
weaker  player  firing  before  the  stronger  one  for  normalized  times  near  t^  =  -1, 
but  also  for  the  weaker  player  firing  after  the  stronger  one  for  large  t. 

Properties  (ixb)  and  (ixc)  combine  to  show  that  the  mode  of  f^  occurs 
at  tg  -  -1  while  that  of  f ^  occurs  at  a  value  of  £  greater  than  that  0 ^ 
which  the  two  players’  hazard  rates  are  equal,  which  is  in  turn  greater  than 
-1  . 

Properties  (ixd)  and  (ixe)  state  that  the  expectation  of  is  finite 

while  that  of  <J>^  is  infinite. 
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Finally,  information  about  the  constants  c.  .  and  the  initial  values 

1  *  J 

i  <  j  can  be  obtained  by  means  of  a  complicated  recursive 
process.  First,  the  above-mentioned  fact  that  v  (-1)-2(J)9  (-1)  =  0  yields, 
from  the  formulas  for  a,  6^  62»  Yp  and  y2: 


v  (-1)  = _ ±11 _ 

m,n'  (m+n-1) 


_  2  (m+n-1)  cm,n  ,  , 

*  Cm,n  +  IStaT  •  ^~T~  *  Vm-l,n^-1^  lf  1 

m-l,n 


<  m  <  n-1 


and  v,  n(-l)  =  2c/  . 

1 >n  1  ,n 

Then,  this  relation  coupled  with  the  fact  from  [3]  and  [6]  that 

vm # i/ c ^  “  2(^2  (t)  >  0  for  t  >  -1  implies,  by  property  (viii) >  that 


^(-1)  -  4>2(-l)  =  zm>n/(81-S2)/2a, 

where  z  is  the  zero  of  the  Weber  parabolic  cylinder  function  D  , ,  with 

IU  »  U  p+1 

p  =  -m/ (m+n) .  Solving  this  relation  for  c  ,  using  the  fact  that  c  >  0. 

m»n  °  'm,n 

Rives 


c  =  -z  A  +  B 
m,n  m,n  m,n  m,n 


where 


Am,n  =  [(n-m) (m+n-1) /2 (m+n)3(m+n-2) ] 1 /2 

E  _  (m+n-1)2  .  Vm,n-1^-1^  vm-l,n^-1^  , 

m’n "  1  ~ — : - — 1  for  1  <  m  <  "-1 

2 (m+n)  m,n-l  m-l,n 


while  for  n  >  2, 


vl,n-l("1) 


“m-1 ,n 


l,n  2 (n+1)  cx  n_x  (n+1)  cl,n-l 


B  -  (n-1)  #  Vn-2,n^“1^ 

n-1 ,n  2n-l  #  c  „ 

n-2,n 


and  lastly  that 


0. 


Thus,  we  have  c 


/ 1/27  at  the  beginning  of  the  recursive  chain. 


1,2 =  2i  ,2 

Next,  we  find  that  c,  »  -z,  A,  +  B, 

l,n  l,n  l,n  l,n 

expresses  c1  in  terms  of  c,  ,  for  n  >  2,  Hence,  the  formula 
i,n  l»n-l 

2 

v,  (-1)  =  2c,  determines  these  initial  values  in  a  simple  recursive  way. 
1 ,n x  l,n 

It  is  clear  then  that  the  values  c  and  v  (-1)  can  eventually  be 

m,n  m,n  7  J 

calculated  in  terms  of  m,n  and  the  zeroes  z  ,  but  simple  formulas  for 

m,n 

those  values  are  not  apparent. 

Thus,  it  is  indeed  fortunate  thav  the  analysis  presented  here  that  is 
germane  to  the  noisy  duel  problem  doe,  not  depend  on  specific  information 
about  the  coefficients  in  (1)  beyond  the  hypotheses  of  the  Theorem.  That 
lack  of  information  is  compensated  for  by  the  condition  that  properties  (i) 
and  (ii)  are  to  hold. 


L 
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Appendix.  The  proof  of  Lemma  7  depends  on  the  following  results,  which  can 
be  derived  by  straightforward  modifications  of  the  argument  given  in  section  9.2 
of  Hille  [2] . 

THEOREM.  Suppose  r  and  p  are  continuous  functions  such  that  r(s)  >  0 
and  p (s)  ^  0  £25  s  real-  Let  yQ  and  be  the  solutions  of 


(Al)  (ry')'  -  py  =  0 

that  satisfy  the  initial  conditions 

yQ(0)  =  0,  Yq (0)  =  1 

and  Yl(0)  =  1,  yj(0)  =  0. 


Then:  (a)  The  limits 


0  =  lim 


S-*x> 


yj  (s) 


=  lim 


yj(s) 
s-xo  y^(s) 


exist ,  and  0  £  y. 

(b)  The  solutions  of  (Al)  passing  through  the  point  (0,1)  that  have, 
y ( s)  >0  and  yf  (s)  £  0  for  s  >  0  are  precisely  those  solutions  y  *  Aypyj 
that  have  0  £  X  £  M-  Moreover ,  every  such  solution  satisfies : 

(b^)  y(s)  >  0  and  yf(s)  £  0  for  all  real 

(bz)  y1 (s)  <  0  over  any  interval  on  which  p(s)  does  not  vani sh 
identically. 

(c)  0  *  y  exactly  when 

Jo  ir(yo)2  +  p(yo)21  =  °°* 

a  sufficient  condition  for  which  is  the  divergence  of 


no 

r 
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20.  ABSTRACT  (Cantinua  an  raeaeaa  *14*  U  woe aaaacy  dnd  Idantlfy  by  block  number) 

It  is  proved  that  a  Riccati  differential  equation  of  a  particular  form  has  a 
unique  solution  satisfying  the  conditions  that  it  is  to  exist  for  Large  values 
of  the  independent  variable  £  and  to  have  its  graph  stay  above  a  certain  line 
for  large  _t.  It  is  then  proved  that  the  solution  exists  for  all  t.  Two  forms 
of  the  solution  are  developed  in  terms  of  the  confluent  hypergeometric  functions. 
An  application  of  these  results  is  made  to  an  asymptotic  stochastic  analysis  of 
a  noisy  duel  problem. 
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